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ABSTRACT 

We present high-resolution X-ray spectra of cloud-shock interaction regions in the eastern and 
northern rims of the Galactic supernova remnant Puppis A, using the Reflection Grating Spectrometer 
onboard the XMM-Newton satellite. A number of emission lines including Ka triplets of He-like N, 
O, and Ne are clearly resolved for the first time. Intensity ratios of forbidden to resonance lines 
in the triplets are found to be higher than predictions by thermal emission models having plausible 
plasma parameters. The anomalous line ratios cannot be reproduced by effects of resonance scattering, 
recombination, or inner-shell ionization processes, but could be explained by charge-exchange emission 
that should arise at interfaces between the cold/warm clouds and the hot plasma. Our observations 
thus provide observational support for charge- exchange X-ray emission in supernova remnants. 
Subject headings: atomic processes — ISM: abundances — ISM: individual objects: Puppis A — ISM: 
supernova remnants — X-rays: ISM 



1. INTRODUCTION 

It has long been thought that the X-ray emis- 
sion from supernova remnants (SNRs) is partly due 
to charge-exchange (CX) pro cesses between neutral s 
and highly-ionized ion s (e.g., iSerlemitsos et al.l Il973f) . 
Wise & Sarazin (1989) first performed a detailed theo- 
retical calculation of CX-induced X-ray emission from 
SNRs. They found that the contribution of CX emission 
is typically 10 -3 to 10 ~ 5 compared with that of ther- 
mal emission (i.e., electron collisional excitation lines). 
This means that CX emission is in general very minor. 
However, as the authors claimed, CX emission could play 
an enhanced role where neutrals are mixed with shocked 
hot gas because of hydrodynamic instabilities in SNRs 
and/or in homogeneities of the interstellar medium. Sub- 
sequently, iLallementl (|2004f ) examined projected emission 
profiles for both CX and thermal emission in SNRs, and 
noted that CX X-ray emission could be comparable with 
thermal emission in thin layers at the SNR edge. It is 
also pointed out that the relative importance of CX X-ray 
emission to thermal emission is proportional to a quan- 
tity, n c V r n~ 2 , with n c being the cloud density, V r the 
relative velocity between neutrals and ions, and n c the 
electron density of the hot plasma. Thus, the higher the 
density contrast (n c /n e ), the stronger the presence of CX 
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X-ray emission would become. 

Observational evidence of CX emission from SNRs was 
found at optical wavelength s more than 30yrs ago (e.g., 
iKirshner fe Chevalier |[l978| ). but is still lacking in the 
X-ray domain. So far, only marginal detections of CX 
X-ray emission were reported in the SMC SNR 1E0102.2- 
7219 (iRasmussen et al.ll2001l) and t he Galactic SNR, the 
Cygnus Loop (|Katsuda et all 1201 1|) . The line intensity 
ratios of (n > 2 — > n = l)/(n = 2 — s> n = 1) in 
H-like O (1E0102.2-7219) or He-like O (Cygnus Loop) 
seem to be higher than expectations from thermal emis- 
sion models, and were interpreted as signatures of CX. 
However, the CX interpretation is still a matter of de- 
bate, because Chandra grati ng observations of 1E0102.2- 
7219 (jFlanagan et al.ll2004j ) did not confirm the XMM- 
Newton result, and uncertainties of Fe L line emissivities 
might cause the spe ctral anomaly in the Cygnus Loop 
(jKatsuda et al.ll20lll ). 

One SNR that might show CX X-ray emission is Pup- 
pis A, a large ( ~50' in diameter) middle-aged (~3700yr: 
I Winkler et ail 1 19881 ) SNR in our Galaxy. Its surface 
brightness is the highest among all Galactic SNRs in 
soft X-rays (E <lkeV). The X-ray image is inhomo- 
geneous with particularly bright eastern and northern 
knots (hereafter, BEK and NK). These regions un- 
doubtedly represent strong interaction s between th e SNR 
blast wave and dense clouds (e.g.. | Petre et al.l 119821 : 
iDubner fe Arn al 1988; Hwang et al.ll2005l) . and thus they 
exhibit ideal conditi ons in which CX X-ray emission 
could be significant (Lall ementl 120041 ). While indica- 
tions of CX signatures have not yet been reported in 
this SNR, such signatures can be best obtained through 
high-resolution X-ray spectroscopy. 

Here, we present high-resolution X-ray spectra of 
the BEK/NK features using XMM-Newton's Reflectio n 
Grating Spectrometer (RGS: Iden Herder et al.l 120011 ). 
From the RGS spectra, forbidden-to-resonance line ratios 
in Hea triplets are found to be anomalously enhanced, 
especially at the BEK. We show that this anomaly could 
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Fig. 1. — Left: Merged XMM-Newton and Chandra image in the energy band 0.5-5 keV. The intensity scales as the square root of the 
surface brightness. The RGS spectral extraction regions are overlaid in white with the region name. The arrows indicate zero points of 
projection profiles shown in Fig. [l] right. Right: X-ray emission profiles along the RGS dispersion axis for the five spectral extraction 
regions. The zero points (x=0) correspond to the direction from which emission is detected at the nominal (no red/blue shift) wavelength 
positions on the RGS detectors. 

generate more accurate RGS response files than the re- 
sponse file for the entire BEK region as we describe be- 
low. We smooth RGS responses originally designed for 
point sources, based on emission profiles along the dis- 
persion direction with the rgsrmf smooth software. For 
the input to the software, we arrange the X-ray images 
such that Chandra covers regions around the BEK/NK 
while XMM-Newton covers the remaining regions, and 
we take account of vignetting effects of XMM-Newton 1 & 
X-ray telescope. We generate energy-dependent RGS 
responses using energy-band images in 0.35-0.4 keV, 
0.4-0.5 keV, 0.5-0.6 keV, 0.6-0.7 keV, 0.7-0.85 keV, 0.85- 
0.97keV, 0.97-f.2keV, and 1.2-1.5 keV. Figure Q] right 
shows the X-ray emission profiles in 0.5-0.6 keV for the 
five regions. These profiles, smoothed by the point 
spread function of the telescope, become the RGS re- 
sponses themselves for the corresponding energy band. 
We note from this figure that the contribution from the 
NK in BEK spectra is almost negligible due to vignetting 
effects (the NK is slightly seen around x=33' in the upper 
four panels in Fig. [1] right). 

For RGS background (BG), we use a blank sky ob- 
servation (Lockman hole: Obs.ID 0147511601) for which 
the RGS spectrum is free from apparent line emission. 
In addition, we need to subtract local BG, since the Vela 
SNR is superposed on Puppis A. To estimate local BG, 
we simulate RGS spectra based on plasma parameters de- 
rived by modeling EPIC spectra outside Puppis A. The 
simulated local BG is added to the blank sky spectrum, 
resulting in a total BG. The total BG is dominated by 
local BG, and is less than 1% of the source emission. As 
for the EPIC BG, we use blank-sky data prepared by 
[Read fc Ponmanl (|2003l ). Since X-ray emission from Vela 
is much weaker than the BEK/NK features, especially in 
the higher energy band, we do not take account of local 
BG for the EPIC spectra. 

Figure [5] shows BG-subtracted spectra for the five re- 
gions. To make the best use of XMM-Newton^ capa- 
bility, three different data sets are analyzed together: 
the first order RGS spectra below 0.6 keV, the second 
order RGS spectra in the 0.65-1.5 keV band, and the 
MOS spectra above 1.1 keV. Note that the second or- 



be naturally caused by the presence of CX emission. 

2. OBSERVATIONS AND SPECTRAL ANALYSIS 

The NK and BEK were observed by XMM-Newton on 
2003 April 17 (Obs.ID 0150150101) and 2003 May 21 
(Obs.IDs 0150150201 and 0150150301), respectively, in 
order to obtain high- resolution spectra using the RGS. 
The dispersion directions of the RGS are 76°. 4 (NK) and 
132°. 3 (BEK) counterclockwise from the north, as shown 
in Fig. [TJ which is an X-ray image of Puppis A gener- 
ated from existing XMM-Newton and Chandra data. In 
this paper, we mainly focus on the RGS data, while we 
also utilize data taken by the European Photon Imaging 
Camera (EPIC: ITurner et al.l[200lt IStruder et afll200l 
to support our RGS analyses. The exposure times after 
removing intervals affected by soft protons are 13.5 ks 
and 20.8 ks for the NK and BEK, respectively. All the 
raw data are processed using version 11.0.0 of the XMM 
Science Analysis Software and the latest calibration data 
files available at the time of the analysis. 

Since the RGS is a slitless spectrometer, off-axis 
sources along the dispersion direction are detected at 
wavelength positions shifted with respect to the on-axis 
source. Spatial displacement of V corresponds to a wave- 
length shift of 0.138 A(or 4eV/arcmin at 0.6 keV) for the 
first spectral order. Because the BEK and the NK are 
fairly compact (2'-3') features with much higher surface 
brightness than their surroundings, the RGS is capable 
of producing high-resolution spectra for them, with an 
order-of-magnitude better resolution than nondispersive 
CCDs. In fact, there are a number of successful RGS 
observations not only of moderately extende d SNRs in 
the LMC/SMC (e.g., iRasmussen et all 120011 ), but also 
of a bright knot a long the northwestern rim of SN 1006 
(|Vink et al.H200l . 

As shown in Fig. [TJ we divide the eastern RGS field 
of view into four sectors spaced by 0'.8 along the cross- 
dispersion axis, while we extract one spectrum from the 
l'.6-width region for the north. In this way, we mit- 
igate the photon number difference between the BEK 
and the NK. Also, slicing the BEK region allows us to 
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der RGS spectra, which show the best spectral resolu- 
tion among all the X-ray spectrometers onboard XMM- 
Newton, do not cover the energy range below 0.65 keV, 
and that above ~lkeV the signal-to-noise ratio of the 
EPIC is much better than that of the RGS. We sum 
spectra from two co-aligned RGS spectrometers, RGS1 
and RGS2, to improve photon statistics. We extract 
MOS spectra from the same sliced regions as the RGS, 
limiting around the BEK or the NK. As expected from 
previous X-ray observations, the RGS spectra are dom- 
inated by a number of emission lines, including C Lya 
@0.375keV, N Hea triplet @0.42keV, N Lya+N He,3 
@0.5keV, O Hea triplet @0.57keV, O Lya @0.654keV, 
O He/3 @0.666keV, Fe L (3s^2p) complex @0.73keV, 
O Ly/3 @0.775keV, Fe L (3d^2p) complex @0.83keV, 
Ne Hea triplet @0.92keV, Mg Hea triplet @1.35keV, 
and Si Hea triplet @1.85keV. In particular, as shown in 
the right (and left) panels of Fig. [21 forbidden and res- 
onance lines in the Hea triplet of O (N and Ne as well) 
are clearly resolved for the first time for the two knots. 
We also note that this is the first clear detection of C 
and N lines from Puppis A. 

We fit the spectra with an absorbed bremsstrahlung 
component plus 49 Gaussian components. Free parame- 
ters are the H column density A^j, the electron temper- 
ature kT C7 and the normalization of the bremsstrahlung 
component. We also allow normalizations of most Gaus- 
sians to vary freely. Since O Ly7 and O Lyi5 lines fall 
in the Fe L line forest, we fix the normalization ratios 
of O Ly 7 /0 Ly/3 and O Lyd/ O Ly/3 to 0.3 and 0.1 as 
are expected at kT c ~0.5keV (|Hwang et all [2005). We 
note that the kT c values of the bremsstrahlung compo- 
nent are mainly determined by continuum emission above 
1.5 keV. Thus, fixing ratios of Gaussian normalizations 
does not affect the kT e measurements. Line centers of 
the most prominent 21 Gaussians are left as free param- 
eters, while those for C He/3 and two Fe L lines (G0 
and D) are systematically shifted with respect to those 
of their neighboring lines (H and E). Other weak Gaus- 
sian s' centers are fixed to the theoretically expected val- 
ues ([Smith et al.ll200Th . The widths of all Gaussians are 
fixed to zero, since significant broadening is not required 
from a statistical point of view; note that apparent line 
widths seen in the RGS spectra in Fig. [5] are due to spa- 
tial extent of the source as described in the second para- 
graph in this section. With this fitting strategy, we ob- 
tain fairly good fits for all the spectra as shown in Fig. O 
The fit results are summarized in Table [TJ where we omit 
weak Gaussians. 

The values o f 7Vh are consistent with recent X-ray 
measurements ([Hwang et all 120051 120011 iKatsuda et all 
120101 ). The kT e values are somewhat lower than previ- 
ous results. However, prominent line intensities, which 
are essential for the discussion below, are not affected 
by the temperature difference, as the fraction of under- 
lying continuum is very small. Line intensity ratios are 
basically reproduced by thermal emission models hav- 
ing plausible plasma parameters (i.e., electron tempera- 
ture of 0.5 keV an d ionization timescale of 10 11 cm~ 3 s: 
iHwang et ai]|200H) . 

On the other hand, as summarized in Table the 

8 We use the s ame labeling for Fe XVII L lines as in 
IGillaspy et"aT] SWHt ). 



intensity ratios of forbidden line to resonance line in 
Hea triplets (hereafter, F/R ratios) are inconsistent 
with model predictions, for which we assume a plane- 
parallel shock model with a distribution of ioniza- 
tion timcscales appropriate for a plane-parallel shock 
(vpshock model in conjunction w ith augmented NEI ver- 
sion 2.0: iBorkowski et all 120011) ; this emission model is 
often used to describe X-ray spectra from Puppis A. 
In calculating the model F/R ratio, we take a reason- 
able range of plasma conditions in the BEK (kT e = 
0.3-0 . 7keV and n c t = 10 10 -5x 10 11 cm' 3 s: IHwang et all 
1200511 given that the RGS spectra are integrated along 
the dispersion direction as well as the line of sight. From 
Tabled we see that the F/R-ratio anomaly for O and 
Ne ions is particularly evident in the E2-E4 regions com- 
pared with that in the El and NK regions, while the 
F/R ratios for N ions are marginally consistent among 
the five regions. The discrepancy can not be solved by 
multiple plasma components, since the model predictions 
already consider various (reasonable) plasma conditions. 
Although the ratios could be reproduced at a very low 
n c t value of ^10 9 cm~ 3 s, in such a low ionization con- 
dition, no Ne-like Fe L lines nor H-like Ne lines can be 
emitted, which is in stark contrast to the observed X-ray 
spectra. Thus, the F/R ratios cannot be reproduced by 
thermal emission models. 

We next investigate whether the anomalous F/R ratio 
is due to reduction of the resonance line or increase of 
the forbidden line, by comparing the observed line in- 
tensities with model expectations. To this end, we fit 
the X-ray spectra (in the E2 region) with an absorbed 
vpshock model, excluding the particularly anomalous 
O Hea triplets. After fitting, we recover the O Hea 
triplet lines to compare the data with the model. We 
find that the intensity of the forbidden line inferred from 
the vpshock model is weaker than the data and that 
the model intensity of the resonance line is stronger. 
This result suggests both reduction of the resonance line 
and increase of the forbidden line. However, we need 
to be careful about limitations of our spectral model- 
ing. This is because the X-ray emission other than the 
O Hea triplets may also not be interpreted in the frame 
of pure thermal emission models, due to contamination 
of CX emission, effects of resonance line scattering, or 
some other processes. If this is the case, the plasma pa- 
rameters including O abundance, the electron tempera- 
ture, and the ionization timescale inferred by our thermal 
emission modeling should have some systematic uncer- 
tainties, leading to an incorrect model intensity of the O 
Hea triplet. Therefore, it is difficult to infer the model 
O Hea intensity accurately enough at this point. We 
will revisit this issue in our future work that will include 
developments of more sophisticated emission models. 

3. DISCUSSIONS 

We have presented high-resolution X-ray spectra of 
the BEK/NK features in the Galactic SNR Puppis A. 
The forbidden and resonance lines in the Hea triplets 
are clearly resolved, and their intensity ratios (F/R) are 
found to be generally higher than predictions from ther- 
mal plasma models. This anomaly is particularly evident 
in the BEK, while the RGS spectrum of the NK as well 
as th e Einstein FPCS spe ctra of the northeastern por- 
tion ([Winkler et al.|[l981allbl ) are closer to predictions of 



4 



Katsuda et al. 



rj 

0.1 
4 
2 

K 



% 
2 

« 
-4 



2 

« 



H 

x 

-2 
-4 



El 


\0A 




wiWw WWW 




0.5 


i 

Energy (keV) 




^^^^^ 








0.5 


1 

Energy (keV) 




E3 

/ 






§ppi| 






0.5 


i 

Energy (keV) 




E4 










0.5 


i 

Energy (keV) 






0.56 0.57 0.58 

Energy (keV) 




0.56 0.57 0.58 

Energy (keV) 



2 

X 



0.56 0.57 0.58 

Energy (keV) 



Energy (keV) 

Fig. 2. — Left: XMM-Newton spectra of the five regions in Fig.[T] The data covering below 0.6 keV, 0.65-1.5 keV, and above 1.1 keV are the 
first order RGS1+2, the second order RGS1+2, and MOS1+2, respectively. These data are simultaneously fitted with a phenomenological 
model (see, text), and the best-fit models are shown in red. Lower panels show the residuals. Right: Close-up RGS spectra for O Heo 
triplets with individual best-fit model components. Forbidden and resonance lines are indicated as "f and "r", respectively. 
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thermal models. 

There are two ways to enhance the F/R ratio: (1) an 
increase of the forbidden line flux or (2) a reduction of 
the resonance line flux. Possible mechanisms for the for- 
mer are either H-like— >-He-like recombination/CX pro- 
cesses or inner-shell ionization of Li-like ions due to the 
higher statistical weight. As for the latter, we consider 
resonance-scattering effects; only resonance lines are 
scattered out of the line of sight because of their higher 
optical depths. Among several SNRs for which high- 
resolution spectroscopy has been performed, anomalous 
F/R ratios have only been reported in tw o LMC SNRs, 
DEM L71 (Ivan der Hevden et aT~l l200l and 0506-68 
a.k.a N23 ([Broersen et al.ll2011f) . In these cases, the au- 
thors attributed the anomaly to recombination and/or 
resonance scattering, paying little or no attention to CX. 
Below we discuss which mechanism mainly works for the 
case of Puppis A. 



Effects of resonance scattering can be evaluated for 
each emission line fr om its optical depth, r (e.g., 
iKaastra fe Mewe1ll995f ). We calculate the value of r us- 
ing the BEK's typical plasma parameters; i.e., ^0.5 so- 
lar abundances, /cT c ~0.5 keV, n^^lO 11 cm~ 3 s, and the 
hydrogen column density of the X-r ay-emitting plasma 
~ 2.4xl0 19 cm~ 2 from n c ~4cm- 3 (jArendt et alJl2010f) 
multiplied by an assumed plasma depth of 3.2 pc, which 
corresponds to the 3' size of t he BEK at a distance 
of 2.2 kpc ([Revnoso et al.l [20031 we refer to this dis- 
tance hereafter). To derive oscillator strengths and 
ion fractions, we utilize the data bases from CHIANTI 
(iDere et al.ll20q9h an d SPEX (specifically, the neij code: 
IKaastra. et al.l I1996T ) . respectively. Then, assuming no 
microturbulence, we calculate r and the escape proba- 
bility, p. The escape probabilities for resonance lines in 
Hca triplets of N, O, and Ne are obtained to be 1, 0.92, 
and 0.92, respectively, while those for forbidden lines are 
all unity because of their very small oscillator strengths. 
Therefore, we find that effects of resonance line scatter- 
ing likely play a modest role, if any, in explaining the 
observed F/R ratios for N and O in all the five regions, 
while the F/R anomaly for Ne could be explained by the 
effects. Whereas optical depths could differ with chang- 
ing the assumed plasma depth and microturbulence ve- 
locity, the difference seems small in our case. For ex- 
ample, in order to reproduce the observed F/R ratio in 
the E3 region, we need a plasma depth of 77 pc (or 2°) 
which is obviously too long for this region. Also, a micro- 
turbulence velocity of 100 km s -1 (which is statistically 
allowed by our RGS spectral analysis) would reduce r to 
^0.12 and increase p to ^0.95 for a reasonable plasma 
depth of 3.2 pc. It should be noted that the Fe L line ra- 
tio, (F+G+H)/C, should be more sensitive to resonance 
scattering than the O VII F/R ratio: we expect t^0.02 
and p~0.99 for F+G+H and r-0.43 and p-0.84 for C 
for the plasma conditions in the BEK. We find in Table [2] 
that the observed (F+G+H)/C ratios are indeed slightly 
higher than those of the vpshock prediction. While the 
measured line ratios suggest the presence of resonance 



scattering for the Fe L lines, the degree of the effect is 
not as significant as what is expected from the anomaly 
of O VII F/R ratios. Therefore, we are required to ex- 
amine other possibilities. 

X-ray-emitting recombining plasm as have been re- 
cently found in several SNRs (e.g., lYamaguchi et al.1 
120091: lUchida et all 120121 and references therein). The 
recombining plasma is characterized by strong radiative 
recombination continua (RRC) and enhanced Lya/Hea 
ratios. We find that signatures of RRCs (i.e., strong re- 
combination edges) are not evident in the RGS and MOS 
spectra in Fig. [2] (and also Suzaku XIS spectra which are 
not shown in this paper), although the Fe L and/or other 
line emission might make it difficult for us to detect such 
signatures. Also, the Lya/Hea ratios for N, O, and Ne 
in the five regions are all well within thermal predic- 
tions as can be seen in Table [2j while the BEK regions 
show somewhat larger ratios than those in the NK. In 
addition, recombination processes would e nhance the Fe 
L line ratio (F+G+H)/(C+D+E) to >25 (jLiedahl et al.l 
1990), which is inconsistent with our RGS measurements 
of 1.2-2.4 in the BEK and the NK. Furthermore, our 
spectral fitting of the combined RGS and MOS spec- 
tra with a recombining plasm a model (the cie model 
in SPEX: IKaastra. et al J [19961) failed to reproduce the 
entire X-ray spectrum; the model requires too low an 
electron temperature to explain emission above 1 keV. 
These investigations indicate that the plasma here is not 
recombining. 

Signatures of inner-shell ionization processes of Li- 
like ions can be found as Li-like satellite lines as well 
as an enhanced forbidden-to-intercombination (F/I) ra- 
tio compared with collisionally excited emission (e.g., 
iPorquet et al.l [20Tot ) . In the RGS spectra, there is no 
indication of satellite lines, however. Also, the measured 
F /I ratios of Hea triplets in Table [2] are all marginally 
consistent with thermal expectations; we need to take 
account of considerable systematic uncertainties on the 
intensity of the intercombination line due to its weakness 
compared to the surrounding forbidden and resonance 
lines. These facts led us to conclude that inner-shell ion- 
ization processes are not working efficiently in this region. 
We note that, whereas line emission from Li-l ike O is seen 
in th e far-ultraviolet spectrum of the BEK ([Blair, et al.l 
1995), the abundance of Li-like O would be small in the 
X-ray-emitting region. 

We next assess the feasibility of the CX s cenario by 
calcul ating the expected CX flux, following lLallemeirn 
(2004). The volume emissivity of CX is expressed as 
Pcx — ocx ni V T . We let ocx be the CX cross 
section between neutral H and ions of interest, njj the 
neutral H density, rii the ion density, and V T the relative 
H-ion velocity We first consider the O VII forbidden 
line, so that acx is 3.3x 10~ 15 cm -2 ([Bodewits. et al.1 
2007$. The value of tt-h can be taken from the den- 
sity of cold/ warm clouds immersed in X-ray-emitting 
plasma, which is ~50cm" 3 (jTeske fe Petrel Il987t ). The 
value of rii is the density of O VIII in the hot plasma , 
which is ~2xlQ- 4 c m~ 3 for the BEK ([Hwang et al.ll2008t 
lArendt et all 120101 ) . We assume that V r = 500kms _1 , 
roughly the shock velocity that can produce the X-ray- 
emitting plasma in the BEK. These parameters give 
P C x^l-7xlO~ 9 photonscm _3 s _1 . The CX-emitting vol- 
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TABLE 1 
Spectral-fit parameters 



Component 


Parameter 






Region 










El 


E2 


E3 


E4 


N 


T7 7; 

Absorption 




9 70-1-0 01 


9 RO-I-O 01 


9 RO-I-O 01 


9 rc;-(-o 01 


9 ccfi-|-0 01 


Bremsstrahlung 


fcT e (keV) 


0.32±0.01 


0.31±0.01 


0.33±0.01 


0.32±0.01 


0.30±0.01 


Normalization" 


5264±93 


9124±133 


6744±117 


7808±117 


5397±90 


Gaussian: C Lya+Si L+S L 


Center (eV) 


366. 8 + o \ 


367.2+Jj'^ 


366.71q3 


365. 4+0 3 


367 




MriT'TTi a 1 1 <7 n 1" 1 <~m 


022+80 


14XB+HX 

iwyj no 




1 79^-1-1 iq 




N Hea (f) 


Center (eV) 


419 2 + S'? 
^ — 0.5 


419.2±0.3 


418.4±0.5 


419.7±0.5 


420 5 +0 ,'? 

^^'"•'-'—0.6 




Normalization" 


100±18 


201±26 


145±26 


243+33 


82+16 


N Hea (i) 


Center (eV) 


426 6 


426 6 


426 6 


426 6 


426 6 


Normalization" 


<15 


57±23 


72±24 


<57 


18+15 


N Hea (r) 


Center (eV) 


429.71q2 


430.ll°| 


429.9±0.5 


429.7±0.5 


430.6+0.5 




Normalization" 


116±17 


223±24 


169±23 


199+27 


72±15 


C Ly/3 


Center (eV) 


434.9+IH 

—0.2 


435.1+°« 

—0.3 


435.2+0-7 

—0.6 


434.7+°I 

— 0.6 


435. 2+°'^ 

—0.6 




Normalization" 


60±14 


87±18 


55±18 


87±22 


51±13 


C Ly 7 


Center (eV) 


459 6 


459 6 


459 6 


459 6 


459 b 


Normalization" 


25±8 


31±10 


24±10 


14±12 


< 8 


N Lya+He^ 


Center (eV) 


499.4±0.2 


500.21° 2 


499.8lgj 


499.6l°2 


499.6l°i 




Normalization" 


143±9 


254±13 


222±12 


252+14 


162±li 


O Hea (f) 


Center (eV) 


559.8l°* 


560.7±0.1 


560.5±0.5 


560.2+0.5 


561.2±0.2 




Normalization" 


705±27 


1167±36 


1042±36 


1119±38 


417±23 


O Hea (i) 


Center (eV) 


566.8lJ° 


568.31°/? 


570.71°-* 


568.51°-* 


567.51?* 




Normalization" 


107±24 


190±30 


298±3i 


255+34 


168±22 


O Hea (r) 


Center (eV) 


573.1±0.2 


574.1±°;1 


573.4±0.1 


572.9±0.1 


574.3!°;! 




Normalization" 


601±25 


771±30 


527±31 


646+33 


492±21 


O Lya 


Center (eV) 


652.7tgi 


653.4±0.2 


653.0±0.2 


652.8±0.2 


653.4±0.4 




Normalization" 


979±38 


1223±47 


1154±46 


1329+51 


631±36 


O He,3 


Center (eV) 


664.4±0.8 


665.0I07 


665.41}'^ 


666.0II4 


667.0±1 




Normalization" 


161±24 


285±29 


197±27 


136+29 


110±17 


O He7 


Center (eV) 


699.0+gg 


700.7±0.6 


698.9l ;i 


698.21J 1 


700.7±1.1 




Normalization" 


68±8' 


74±10 


65±11 


74±13 


42±9 


Fe L (G+H) 


Center c (eV) 


724.9l°* 


725.ll°i 


724.9l°;i 


725.2l°; 3 


725.2l° 4 




Normalization" 


2471JI 


362l^ 


33lUl 


332+j° 


156±i| 


Fe L (F) 


Center (eV) 


739. 3+n? 

— U.O 


738.5l°'2, 


739.3+°'® 

— U.o 


738.1+°® 

— U.o 


738.9 + n ? 

— U.o 




Normalization" 


89±11 


156±13 


94±12 


99±12 


65±8 


O Ly/3 


Center (eV) 


773.5+°,;? 


773.8±0.3 


773.31°-* 


772.91°-* 


773.7±0.5 




Normalization" 


151±7 


187±8 


160±7 


169+8 


97±6 


Fe L (D+E) 


Center c (eV) 


810.6±0.8 


810.1±0.2 


810.0±0.4 


810.0+0.4 


810.1±0.7 


Normalization" 


74+1° 


no±iS 


98tlo 




50!^ 


Fe L (C) 


Center (eV) 


826.61°/* 


826.1±0.5 


825.6±0.5 


826.1±0.5 


826.2±g| 




Normalization" 


105±6 


173±0.8 


135±7 


152+8 


78±5' 


Ne Hea (f)+Fe L 


Center (eV) 


902.0±0.6 


903.3lo2 


902.61°;* 


902.61°;* 


903.51°° 




Normalization" 


122±6 


185±7 


156±6 


175+7 


75±4 


Ne Hea (i)+Fe L 


Center (eV) 


915 6 


915 6 


915 6 


915 6 


915 6 


Neormalization" 


63±6 


53±7 


68±6 


50+7 


18±4 


Ne Hea (r)+Fe L 


Center (eV) 


921.2±g;| 


921. 6lg 3 


921.2l°| 


921.8l°3 


922.4±0.2 




Normalization" 


141±6 


224±7 


168±6 


188+7 


111±5 


Ne Lya+Fc L 


Center (eV) 


1023.6+J-4 


1022.5 +l j , -2 


1020.5+°/4 


1022.1+5 4 


1020.41J;* 




Normalization" 


83±5 ' 


97±6 ' 


96±5 ' 


102+6' 


35±3 


rVd-o.f. 




1423/968 


1740/1151 


1698/1114 


1664/1181 


949/612 



NOTE. — a In units of 10 4 photons cm -2 s - 1 . 
b Fixod values. 

c Linc centers of H and E are shown, while those of G and D are systematically shifted by +2eV from H and E, respectively. 



time is calculated by the thickness of the CX layer, lex, 
times the interface area between the clouds and the 
hot plasma. The value of lex is equated to the mean 
free path for H-proton CX. Thus, lex is of the order 
of l/(7?T,p, with a being the H-proton CX cross section 
(10~ 15 cm -2 : iMcClurel Il966t ) and n p being the proton 
density in the plasma (4cm~ 3 : lArendt et al.ll2010T ) . Con- 
sidering that only ~30% neutral H can charge- excha nge 
before getting collisionally ionized (|Lallementl 120041) for 
the BEK's plasma condition, we obtain an effective CX 



thickness of ^7.5xl0 13 cm. An assumed cylindrical in- 
terface with a diameter of 3' would give a surface area of 
the CX layer of ^3xl0 37 cm 2 for each region, resulting 
in a CX-emitting volume of ~2.3x 10 51 cm 3 . With P C x 
and the emitting volume estimated, and an assumption 
of a uniform PcX: w e obtain an unabsorbed CX flux for 
the O VII forbidden line of ~6xl0 -3 photons cm~ 2 s _1 
for each region in the BEK. Next, we repeat the cal- 
culation for forbidden lines of He-like N and Ne. Us- 
ing the relation that acx is proportional to the atomic 
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TABLE 2 
Line intensity ratios 



Element 


Thermal predictions" 






Region 










El 


E2 


E3 


E4 


N 


N F/R 


0.38-0.52 


0.68±0.15 


0.90±0.15 


0.86±0.19 


1.22±0.23 


1.14±0.33 


O F/R 


0.42-0.56 


1.17±0.07 


1.51±0.08 


1.98±0.13 


1.73±0.11 


0.85±0.06 


Ne b F/R 


0.44-0.71 


0.87±0.06 


0.83±0.04 


0.93±0.05 


0.93±0.05 


0.68±0.05 


N F/I 


5.8-6.4 


>7.9 


3.49±1.46 


2.02±0.76 


>3.7 


4.48±3.73 


O F/I 


4.2-5.1 


6.58±1.49 


6.14±0.99 


3.49±0.38 


4.39±0.60 


2.48±0.35 


Ne 6 F/I 


3.3-4.9 


1.92±0.20 


3.45±0.47 


2.30±0.22 


3.46±0.50 


4.01±0.87 


(N Lya + N He/3) / N Hea 


0.12-0.9 


0.58±0.07 


0.53±0.05 


0.58±0.06 


0.54±0.06 


0.94±0.13 


O Lya / O Hea 


0.01-1.03 


0.69±0.03 


0.57±0.03 


0.62±0.03 


0.66±0.03 


0.59±0.04 


Ne 6 Lya / Ne 6 Hea 


0-0.83 


0.25±0.02 


0.21±0.01 


0.24±0.01 


0.25±0.02 


0.17±0.02 


Fe L (F+G+H)/C 


1.4-2.6 


3.19±0.26 


2.98±0.14 


3.15±0.23 


2.84±0.21 


2.83±0.26 


NOTE. — a Predictions by the vpshock 1 


nodcl at (fcT e , net) = (0.3-0.7 kcV, 


10 10 -5xl0 11 cr 


a 3 BlHwane et al. 


11200511. fa The observed v 


aluea could be affect 


cd by Fe L lines. 



number, Z, and the ion number densities of tinvi 
8 x 10~ 6 cm~ 3 and riNoix ~ 4 x 10~ 5 cm" 



based on plasma parameters in iHwang et al.l 120051 ) , we 
estimate the unabsorbed CX fluxes of the N VI and Ne 
IX forbidden lines to be ~ 3x 10 -4 photons cm -2 s _1 and 
~ 2 x 10~ 3 photons cm -2 s _1 , respectively. 

We should keep in mind that the uncertainty of the 
CX flux estimate is fairly large. For example, we im- 
plicitly assume that hydrogen in the cold/warm clouds 
is entirely neutral due to rapid recombination caused by 
radiative cooling in the dense cloud. The hydrogen must 
be partially ionized, however, resulting in a lower CX 
flux. On the other hand, hot neutral hydrogen formed 
after H-proton CX reactions, dust destruction, and un- 
shocked cold neutrals going through the (collisionless) 
shock would all provide additional neutrals in the hot 
plasma, increasing the CX flux. Also, the area of the 
interface region is quite uncertain. These factors suggest 
that our CX flux estimate should only be considered ac- 
curate to an order of magnitude. 

To compare the expected CX flux with our mea- 
surements in Table [TJ we need to correct for the 
spatial-integration factor. This is because the expected 
flux is estimated for the BEK feature, whereas the fluxes 
in Table [T] are integrated along the RGS dispersion direc- 
tion. To estimate the correction (reduction) factors, we 
extract EPIC spectra from 3'-long regions including the 
BEK, and normalize the RGS flux to equalize the EPIC 
flux in each region. In this way, the correction factors 
are derived to be 3.7-4.5 for the four BEK regions. The 
fluxes of the N, O, and Ne Hea forbidden lines are, 
respectively, derived to be ^5 x I0~ 3 photons cm' 2 s _1 , 
~3x 10 -2 photons cm -2 s _1 , and 
~5 x!0~ 3 photons cm~ 2 s _1 for each region in the 
BEK. We thus find order-of-magnitude agreement 
with the expected CX fluxes. Given the considerable 
uncertainty of the expected CX flux, we cannot rule 
out the CX scenario with this level of agreement. In 
addition, there are two pieces of observational evidence 
supporting the CX scenario. First, optically-emitting 
clouds, which we think are the major electron donors, 
are seen only in the BEK, whereas such clouds are not 
detected in the other regions including the NK and 
the Einstein FPCS field of view. Second, the N/O 
and Ne/O abundance ratios based on CX fluxes in the 
Hea triplets for the BEK region are consistent with 
previous X-ray measurements (N/0~1 and Ne/0~2, 



res pectively) . Here , we applied the same method used 
in iLiu et all fj2011h who modeled Hea triplets with a 
thermal emission model plus Gaussians to represent 
the CX emission for the starburst galaxy M82. In this 
context, we conclude that CX emission is a promising 
mechanism for explaining the anomalous F/R ratios 
observed in Puppis A. 

One might expect to find other spectral signatures of 
CX in addition to the Hea line ratios. However, CX 
spectral properties a re strongly dependent on V r an d the 
target neutrals (e.g.. lBei crsdorfcr. et al. 2001, 20031) , and 
investigation of CX emission is still ongoing. Therefore, 
further discussion calls for more sophisticated CX emis- 
sion modeling, which is beyond the scope of this paper 
and is left as future work. 

4. CONCLUSIONS 

High-resolution X-ray spectra of the cloud-shock in- 
teraction regions, the BEK and the NK, in Puppis A 
have revealed anomalous Hea triplet ratios: in particu- 
lar, the O Hea F/R line ratios are found to be ~2 in the 
BEK. This anomalous ratio can be naturally interpreted 
as the result of the presence of CX emission, although 
resonance-scattering effects might be non- negligible. 

Future observations wi th the non-dispersiv e Soft X- 
ray Spectrometer fS XS: iMitsuda et al.l 1201 Of) onboard 
the Astro-H satellite (|Takahashi et al.l |20101T will allow 
for further tests for the CX scenario and other possibil- 
ities. While the RGS can perform high-resolution spec- 
troscopy only for bright knotty features, the SXS with 
expected spectral res olution of 4 eV FWHM and spatial 
resolution of I ' HPD (jSerlemitsos et al.ll2010h will enable 
us to reveal a spatial distribution of the F/R ratios of the 
knotty features and their surroundings. The CX scenario 
can then be examined by seeing if there is a spatial cor- 
relation between the F/R ratios and optical clouds which 
we here consider to be main electron donors. We expect 
that other large Galactic SNRs, especially the Cygnus 
Loop introduced in Section I, would be also promising 
sites to check the presence of CX X-ray emission. 
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